CMU-HEPOO-03 
UTPT-00-07 
August 2000 

NRQCD Analysis of Bottomonium Production at the Tevatron 



Eric Braaten 

Physics Department, Ohio State University, Columbus, OH 43210, USA 



Sean Fleming 

Physics Department, University of Toronto, Ontario, M5S 1A7, Canada 



Adam K. Leibovich 
Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213, USA 



Abstract 

Recent data from the CDF collaboration on the production of spin-triplet bot- 
tomonium states at the Tevatron pp collider are analyzed within the NRQCD 
factorization formalism. The color-singlet matrix elements are determined 
from electromagnetic decays and from potential models. The color-octet ma- 
trix elements are determined by fitting the CDF data on the cross sections 
for T(15'), T(25), and T(35) at large pT and the fractions of T(15') coming 
from Xbi^P) ^-nd Xbi'^P)- We use the resulting matrix elements to predict the 
cross sections at the Tevatron for the spin-singlet states rif,{nS) and hb{nP). 
We argue that r/f,(lS') should be observable in Run II through the decay 
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I. INTRODUCTION 



The NRQCD factorization formalism provides a systematic framework for analyzing the 
inclusive production of heavy quarkonium |jl|] . Long-distance effects involving the binding of 
a heavy quark-antiquark pair into quarkonium are factored into parameters called NRQCD 
matrix elements. These nonperturbative parameters are universal, so values extracted from 
one high energy physics experiment can be used to predict the production rate in others. The 
NRQCD matrix elements scale as definite powers of f , where v is the typical relative velocity 
of the heavy quark. The NRQCD factorization approach becomes phenomenologically useful 
upon truncating the expansion in v so as to reduce the independent NRQCD matrix elements 
to a manageable number. The truncation is most reliable for the heaviest quarkonium states, 
namely the hb system for which is roughly 1/10. 

The most abundant source of data on bottomonium production is the Tevatron pp col- 
lider. In Run lA of the Tevatron, the CDF collaboration was able to resolve the individual 
S-wave bottomonium states T(IS'), T(2S'), and T(3S') and measure their production cross 
sections . An analysis of the CDF data within the NRQCD factorization formalism was 
carried out by Cho and Leibovich [Q. The analysis is complicated by the production of 
P-wave bottomonium states that subsequently make transitions to S-wave states. Cho and 
Leibovich found that the CDF data was insufficient to determine all the important NRQCD 
matrix elements and they had to make educated guesses for some of them. 

The CDF collaboration has recently analyzed the data on bottomonium production from 
Run IB at the Tevatron. In addition to much higher statistics on the cross sections for T(IS'), 
T(2S'), and T(3S') they also have results on the production of the P-wave states ^^(IP) 
and Xfe(2P) 0. The high quality of the new CDF data justifies an updated theoretical 
analysis, with careful attention to the experimental and theoretical errors. 

In this paper, we present a quantitative analysis of the new CDF data on bottomonium 
production within the NRQCD factorization formalism. The color-singlet NRQCD matrix 
elements for S-wave states are determined from their electromagnetic decays, while those 
for P-wave states are estimated from potential models. The color-octet NRQCD matrix 
elements are determined by fitting the CDF data, taking full account of the feeddown from 
transitions between bottomonium states. The resulting values of the matrix elements are 
used to predict the cross sections for the spin-triplet and spin-singlet bottomonium states 
in Run II of the Tevatron. 

II. NRQCD MATRIX ELEMENTS 

The NRQCD factorization approach provides a model-independent framework for ana- 
lyzing the inclusive production of heavy quarkonium []T|. The factorization formula for the 
differential cross section for the inclusive production of a bottomonium state H of momentum 
P has the schematic form 

da[H{P)] = ^rfa[66(n,P)](O^H), (1) 

n 

where the sum extends over both color-singlet and color-octet and over all angular momen- 
tum channels for the hb pair. The hh cross sections, which are independent of the bottomo- 
nium state if, can be calculated using perturbative QCD. All dependence on the state H is 
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factored into parameters {0^{n)) called NRQCD matrix elements. These phenomenological 
parameters can be expressed as matrix elements in an effective field theory called nonrel- 
ativistic QCD (NRQCD). A nonperturbative analysis of NRQCD reveals how the various 
matrix elements scale with the typical relative velocity v of the heavy quark in quarkonium. 
Spin symmetry, which is an approximate symmetry of QCD, also gives relations between 
various matrix elements. 

The relative importance of the terms in the factorization formula (|T]) depends on the 
size of the hh cross sections and on the size of the matrix elements. According to the 
velocity-scaling rules, the most important matrix element for direct T(15') production is the 
color-singlet parameter {Ol^^^^^ i^Si)). The spin-symmetry relations can be used to reduce 
the next most important matrix elements to three color-octet parameters: {(D^'^^^^i^Si))^ 

(Oj^^'^'*(^5'o)), and (Oj^^'^^(^Po))- These color-octet matrix elements are important, because 
the cross sections for producing color-octet hh pairs can be much larger than for color-singlet 
hh pairs. There are analogous matrix elements that describe the direct production of T(25') 
and T(35'). The NRQCD factorization formula (|I|) for direct T{nS) production reduces to 

+da[hbsCSo)]{Oj^'''\'So)) + (Y.i^J+l)da[hbsCPj)]^ {Oj^'''\'Po)) . (2) 

The factor of 2 J + 1 in the last term comes from using a spin-symmetry relation to eliminate 
{Oj^-'\^Pj)) in favor of (Og^^-^^'Po)). 

The most important matrix elements for the direct production of the P-wave states 
X6j(l-P), J = 0, 1, 2, can be reduced to a color-singlet parameter {Oi''°^^^\^Po)) and a single 
color-octet parameter {Og''°^^^\^Si)). There are analogous matrix elements that describe 
the direct production of Xbj{'^P) and Xbj(3-P). The NRQCD factorization formula for 
direct Xbj{nP) production reduces to 

da[xu{nP)] = da[hh,{' P,j)]{Of'^''''\' Pj)) + {2J +l)da[hU'S,)]{Or^''^\'S,)). (3) 

In the last term, the factor of 2 J-(- 1 comes from using a spin-symmetry relation to eliminate 
(Og"*-"'^^(^5'i)) in favor of (Og''°*^""^^(^S'i)). We can also use a spin-symmetry relation to 
replace {Of-'^''^\^Pj)) in the first term by (2 J + l){Of°^''^\^Po)). The matrix elements 
for T{nS) and Xbj{nP) enumerated above should be sufficient for a quantitative description 
of the production of S-wave and P-wave bottomonium states. 

The NRQCD factorization formula gives the cross section for the direct production of 
a given bottomonium state. The cross sections that are most easily measured in experi- 
ments are inclusive cross sections that include contributions from the direct production of 
higher bottomonium states which subsequently decay into the given state. For example, the 
feeddown from Xbi^P), T(25'), and Xbi'^P) accounts for roughly 27%, 11%, and 11% of the 
T(IS') cross section, respectively 0. Taking into account the feeddown from higher T{mS) 
and Xbj{f^P) states, the cross section for inclusive T{nS) production can be written 

dcr[T(n^)]i„, = da[hU'S,)]{0^{^S^))'^}f^ + Y.da[hU^Pj)]{0^{'Pj))'l}f^ 

J 

+da[hh('S,)]{0,{'S,))'l}f'^ + da[hh{'S,)]{0,{'S,))^^f^ 
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TABLE I. Inclusive branching fractions Bh^h' (in %) for transitions between spin-triplet 
bottomonium states. The entries "0 ?" in the first row indicate that the feeddown from X6j(3P) 
is neglected in our analysis. 



+ (^X:(2^ + l)dam'Pj)]^ {OsCPo))I^f\ (4) 

where the "inclusive NRQCD matrix elements" are 

{0[n])l!:''^ = E5^^T(„5)(0^M). (5) 

H 

The sum over H includes T{nS) and all higher bottomonium states that can make transitions 
to T{nS). The coefficient Bh-*h' is the inclusive branching fraction for H to decay into H' . 
By convention, we define Bh~,h = 1- The inclusive branching fraction for the observed 
bottomonium states are collected in Table |. These numbers were obtained by combining 
the measured branching fractions for the exclusive decays T{nS) — > Xfej(mP)-|-7, Xbj{nP) 
T(mS') +7, and T{nS) — > T(mS') + 7r7r, with the exception of -Bt{35)^t(i5), which is a direct 
measurement [^]. 

In Table |, we have not included the spin-singlet states 7]i){nS) and hf){nP), which have 
yet to be observed. Transitions between spin-singlet and spin-triplet states are suppressed, 
because they proceed through magnetic AS* = 1 transitions. The rates for AS = 1 tran- 
sitions are suppressed relative to those for AS* = transitions by a factor of v"^, which is 
roughly an order of magnitude. The branching fractions for r]i,{2S) into other bottomonium 
states are further suppressed by its large annihilation width into two gluons. Quantitative 
estimates of the electromagnetic and hadronic transition rates are given in Refs. They 
support the conclusion that the branching fractions for decays of spin-singlet states into 
spin-triplet states can be neglected. 

The tiny branching fractions in Table | for the transition Xbji'^P) ~* Xbj'i^P) are the 
contributions from the double radiative transitions via T{2S). We have not included the 
contributions from the two-pion decays XbjC^P) XbJ'i^P) + tttt, which have not been 
observed. We can estimate their magnitude by observing that the rates for T{3S) — >■ T(2S')-|- 
717T and T{3S) —>■ T{2S) +77 are equal to within experimental errors. Since the phase space 
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available for the transitions xtjC^P) — > XbJ'i^P) is similar to that for T{3S) T{2S), 
we expect the rate for Xbj{2P) — >• Xbj/(1P) + tttt to be comparable to that for Xbji'^P) 
Xbj'i^P) + 77- Including the effects of two-pion transitions could increase the branching 
fraction for Xbji'^P) Xbj'i^P) by a factor of 2 or 3, but since the values of BxbjC^PJ^Xbj'i^P) 
in Table | are all less than 1.5%, they should still be negligible. 

As indicated by the entries "0 ?" in the first row of Table |I|, we neglect the feeddown 
from the Xb(3-P) states, which have not been observed. A naive extrapolation from the other 
entries of the Table suggest that the branching fractions for Xbii'^P) and Xfe2(3-P) into T{3S) 
could be about 12%, while their branching fractions into T(15') could be about 7%. These 
are small enough that they would not have a significant effect on our analysis. We have also 
neglected the feeddown from D-wave states. 



III. PARTON DIFFERENTIAL CROSS SECTIONS 

In hadron collisions, bottomonium with transverse momentum px of order rrih or larger 
is produced, at leading order in a^, by parton fusion processes ij ^ bb + k. The differential 
cross section for producing a bottomonium state H with momentum P can be expressed in 
the schematic form 

(P)]fosion = fi/p ® fj/p ® da[tj ^ bb{P, n) + k] (0^(n)), (6) 

where there is an implied sum over the partons i,j, k and over the bb channels n. 

The order-a^ fusion cross section in gives a good first approximation only if the 
transverse momentum is not too much larger or too much smaller than rrib- For pt S> rUb, 
the order-a^ fusion cross section for the channel bbsi^Si) has the scaling behavior da /dp\ ~ 
1/pfn, while all other channels are suppressed by powers of m1/p\ at leading order. Par- 
ton processes with scaling behavior are called fragmentation processes. The fragmentation 
contributions to hb channels other than hb^i^Si) enter at higher order in a^. The order-a^ 
fusion cross sections therefore underestimate the bb cross section in these channels at large 
Pt- However, the CDF data on bottomonium production extends only out to px = 20 
GeV m, which is not large enough for fragmentation effects to dominate. In extracting the 
NRQCD matrix elements from that data, it should therefore be sufficient to use the fusion 
cross section (j^). 

The order-a^ fusion cross section in (j^) also fails at small pt- For some bb channels, 
including hb^i^So) and &&8(^-Po,2), there is an order-a^ fusion cross section from the parton 
process ij bb., which produces a hb pair with px = 0. In these channels, the order-a^ 
fusion cross sections da/dp^ diverge like 1/pf. as ^ 0. The divergence in the integral of 
the cross section for ij —>■ bb + k is cancelled by the radiative corrections to the cross section 
for ij — »• bb, so that the cross section integrated over p^ is finite order by order in a^. In 
order to obtain a smooth prediction for da/dp^ in the small px region, it is necessary to 
resum higher order corrections involving soft-gluon radiation. This resummation will have 
a significant effect on the shape of the px distribution, and therefore on the values of the 
NRQCD matrix elements used to fit that distribution. We will avoid the complications 
due to soft-gluon radiation at small px by using only the data from > 8 GeV to fit the 
NRQCD matrix elements. 
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We proceed to describe each of the factors in the fusion cross section (|) in more detaiL 
We include the contributions from the following combinations of colliding partons: ij = 
99i 91i 9^1 QQy where q = u,d,s,c. The parton distributions fi/p{xi,^F) and fj/p{x2, fJ^p), 
which specify the momenta of the colliding partons, depend on a factorization scale //p. 
We will consider the CTEQ5L and MRST98LO parton distribution functions. They are 
both obtained from leading order analyses, and thus can be used consistently with leading 
order parton cross sections. Explicit expressions for the parton differential cross sections 
da are given in Ref. and in Ref. 0. They are proportional to a^{fiji), where /i/j is the 
renormalization scale, and they depend on the mass rrib of the bottom quark. As part of 
the theoretical error, we will allow fip and to vary by factors of 2 from the central values 
fiT = if^b +PtY^'^- This central value interpolates between half the partonic invariant mass 
at bb threshold and half the partonic invariant mass at large pt and central rapidity. 

The cross sections also depend on two fundamental QCD parameters: as and m^. We take 
the QCD couphng constant ^^(/i) to run according to the one-loop formula, with the bound- 
ary value appropriate to the parton distribution function that is being used. For CTEQ5L, 
the coupling constant satisfies as{Mz) = 0.127 and as{mf,) = 0.232. For MRST98LO, the 
coupling constant satisfies as{Mz) = 0.125 and as{mb) = 0.226. 

The other fundamental QCD parameter in our calculation is the bottom quark mass rUb- 
There have been several recent determinations of mb using sum rules calculated to next-to- 
next-to-leading order accuracy with nonrelativistic resummation |10|. A useful summary of 



these results is given by Beneke in Ref. |]TT|. The value of the pole mass is rather unstable 
under radiative corrections compared to short-distance definitions of the mass, such as the 
running mass evaluated at its own scale, rhb = mb{fhb). Beneke's best estimate for this 
mass is fhb = 4.23 ± 0.08 GeV. The definition of rhb is purely mathematical in character 
and not related to any physical thresholds involving the b quark. Two definitions that 
are also relatively stable under radiative corrections and whose definitions are related to 
thresholds in the bottomonium system are the IS mass, which is the perturbative energy of 
the lowest bound state, and the PS mass, which is the sum of the pole mass and some energy 
related to the potential between the b and b. Beneke's best estimates for these masses are 
mb,is = 4.77 ± 0.11 GeV and mb,Psi2GeV) = 4.57 ± 0.10 GeV. The relation between the 
two is given by a power series in a^: 

rrib.psiN = 'n^b,is ^ /i + 0(aj, (7) 

OTT 

The difference between Beneke's values for rub^is and rrib^ps is mostly accounted for by the 
order-tts correction. We will choose the IS mass as our prescription for the b quark mass. 
Beneke's central value for the PS mass differs by 2 standard deviations from the IS mass. 
This difference should not be regarded as an ambiguity in the quark mass, because it could 
not be eliminated by a more precise determination of m^. Instead its effects on the cross 
section could be decreased by calculating the next-to-leading order radiative correction to 
the parton cross sections. The uncertainty due to different prescriptions for the quark mass 
can therefore be regarded as part of the error due to radiative corrections. 
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^qX6o('^-P)|-3 

potential models 


lattice 


T(35) 
T{2S) 
T{1S) 


4.3 ± 0.9 
4.5 ± 0.7 
10.9 ± 1.6 


3.7 ± 1.5 
5.0 ± 1.8 
10.8 ± 5.5 


9.6 ± ? 
3.6 ± ? 
7.6 ± ? 


X6(3P) 
X6(2P) 


2.7 ± 0.7 
2.6 ± 0.5 
2.4 ± 0.4 


1.5 ± ? 



TABLE II. Direct color-singlet matrix elements for T and Xb states ((Oi(^S'i)) in units of 
GeV^ (Oi(^Pj)) in units of GeV^). 

IV. COLOR-SINGLET MATRIX ELEMENTS 

The color-singlet matrix elements for T{nS) can be determined phenomenologically from 
its decay rate into a lepton pair. The electronic decay rate of the T{nS), including the QCD 
radiative correction of order a, and the first relativistic correction of order v'^, is 



T[T{nS) e+e 



9ml 



8 _ 1 Mr(nS) - "^rrib 
3 TT 3 



where Cb = —1/3 is the bottom quark charge and mj, is the IS mass. The relativistic 
correction was first expressed in terms of MY(nS) ~ '^^b by Gremm and Kapustin [^]. The 
vacuum saturation approximation, which is accurate up to corrections of order f ^, has been 
used to express the NRQCD matrix element that enters naturally in annihilation rates in 
terms of the corresponding production matrix element {Oj^"'^\^Si)). It has also been used 
to express the radiative and relativistic correction factor as a square. Setting mb = 4.77 GeV 
and as{m,b) = 0.22 and using the measured value for the decay rates, we obtain the values for 
the color-singlet matrix elements in Table |T[ In addition to the experimental errors in the 
decay rates, there are theoretical errors from relativistic corrections and from perturbative 
corrections. As a measure of the relativistic error of order f ^, we take the square of the largest 
of the order-t>^ corrections for the three T{nS) states: [(Mx(35) — 2m;,)/(3mfc)]^ ^ 0.2%. As 
a measure of the perturbative error from higher orders in a^, we take the square of the order- 
as correction in (|^): [16Q;s/(37r)]^ ^ 14%. The error bars quoted in Table |T|are obtained by 
combining the experimental, relativistic, and perturbative errors in quadrature. The error 
bars are dominated by the 14% perturbative error, except in the case of the T(3S') for which 
the experimental error is 16%. The values for {Oj^''^\^Si)) in Table g are larger by a factor 
of 3 than those given in Table I of the first paper in Ref. [0] because of a normalization error 
in the Table. This normalization error did not appear in the cross sections and therefore 
did not affect the results. 

There is no data that can be used for phenomenological determinations of the color- 
singlet matrix elements for the P-wave states. However the color-singlet matrix elements 
for both the S-wave and P-wave states can be estimated using wavefunctions from potential 
models. Using the vaccuum-saturation approximation, the color-singlet matrix element for 
T{nS) can be expressed in terms of its radial wavefunction at the origin, while that for 
Xbj{nP) can be expressed in terms of the derivative of its radial wavefunction at the origin: 
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T{3S) 


4.3 ± 0.9 


? 


? 


? 


T{2S) 


5.0 ± 0.7 


0.12 ± 0.06 


0.55 ± 0.15 


0.42 ± 0.10 


T(15) 


12.8 ± 1.6 


< 0.2 


1.23 ± 0.25 


0.84 ± 0.15 



TABLE III. Inclusive color-singlet matrix elements for T states ((Oi(^«S'i)) in units of GeV^, 
(Oi(3Pj)) in units of GeV^). 



{01 



T{n5)/3 
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'S'l)) ~ ^I^T{n5) 



(9) 
(10) 



Eichten and Quigg have tabulated the radial wavefunctions and their derivatives at the 
origin for 4 potential models that reproduce the observed bottomonium spectrum [0. As 
estimates of the color-singlet matrix elements (|^) and (p!OD, we take the mean values from 
the 4 potential models in Ref. |T^. The mean values are tabulated in Table The errors 
are the root-mean-square deviations of the 4 potential-model values. The potential model 
estimates of {Oi^'^^\'^Si)) are consistent with the phenomenological values, but have larger 
error bars. This gives us some confidence in the potential-model estimates for {Oi''°^^^\^ Pq)) 
in Table ||. These values are consistent within errors with those used in the analysis of Ref. 



The most accurate determination of the color-singlet matrix elements for the lowest 
bottomonium states may eventually come from lattice gauge theory. The corresponding 
annihilation matrix elements can be readily calculated using lattice simulations of NRQCD 
I^ . The NRQCD collaboration has calculated the wavefunctions for the lowest S-wave 
and P-wave states [|1^] of bottomonium. Inserting these wavefunctions into the 



states Ip^ ^ 

expressions in (^) and (|iy), we obtain the estimates of the color-singlet matrix elements 
in Table [I^. The largest errors in the lattice calculations come from matching of lattice 
NRQCD operators with continuum NRQCD operators and from the omission of dynamical 
light quarks. Both errors could be as large as 25% in the present calculations. It is premature 
to quote error bars for the lattice gauge theory results in Table ||. 

We will adopt the phenomenological values of {Oi^"'^\^Si)) in Tableland the potential- 
model values for (Oj^'°*^"^''(^5'i)). Using the branching fractions in Table |I|, we can form the 
linear combinations that appear in the expressions for the inclusive T[nS) cross sections. 
These are tabulated in Table As indicated by the zeros in the first row, we neglect 
feeddown from the X6(3-P) states. 



V. COLOR-OCTET MATRIX ELEMENTS 

The color-octet NRQCD matrix elements are phenomenological parameters that can only 
be determined from experimental data. We first extract the inclusive color-octet matrix 
elements for T{nS) from the CDF measurements of the inclusive T{nS) cross sections. We 
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then extract direct color-octet matrix elements for Xbj{nP) from the CDF measurements of 
the fraction of T(IS') coming from x^'s. This gives us enough information to determine the 
direct color-octet matrix elements for T{nS) . 



A. Inclusive Matrix Elements for S-waves 



The inclusive T{nS) cross sections depend linearly on the inclusive matrix elements 
defined in (|]). The inclusive color-singlet matrix elements are given in Table |Tl|. We can 
extract the inclusive color-octet matrix elements from the CDF measurements of the T{nS) 
cross sections [|]. The differential cross sections integrated over rapidities \y\ < 0.4 have been 
measured out to pt = 20 GeV @]. The CDF data on Bda/dpx for T{nS), where B is the 
branching fraction for T{nS) f^~^fJ'~ , are shown in Figures 1, 2, and 3 for T(15'), T(25'), 
and T(35'), respectively. We avoid the problem of carrying out soft-gluon resummation to 
determine the shapes of the theoretical pt distributions at low pt by using only the data 
from Pt > 8 GeV to fit the color-octet matrix elements. This leaves 5 pt bins for T(IS') 
and 3 pt bins each for T(2S') and T(3S'). 

The inclusive T{nS) cross sections depend on the inclusive color-octet matrix elements 
through the linear combination [(O8(^5'o)) + m {Os(^Po))/ml] +n {Os{^Si)), where m varies 
from 4.6 at = 8 GeV to 3.4 at pt = 20 GeV, while n varies from 1.0 at pt = 8 GeV 
to 6.3 at Pt = 20 GeV. The parameters (08(^5*0)) and {Os{^Po)) can not be determined 
independently, because the corresponding parton cross sections have similar dependences 
on Pt- We therefore carry out our analysis under the two extreme assumptions that either 
the (08(^>S'o)) term or the {Os{^Po)) term dominates and that the other can be neglected. 
Assuming both matrix elements are positive, the truth will be somewhere in between the two 
extremes. We will take the difference between the two extremes as part of the theoretical 
error. 

For a given choice of the parton distributions and the scales fip and fin, we can integrate 
the T{nS) differential cross section (^) over \y\ < 0.4 and over each pt bin. We determine 
the best fits for the color-octet matrix elements by minimizing the associated with the 
sum over pt bins. Fixing fip = fiji = fiT and assuming that {Os{^Po)) is negligible, we obtain 
the values of {Os{^Si)) and {Osi^So)) in the 1st and 3rd columns of Table Repeating 
the analysis but assuming that {Ogi^So)) is negligible, we obtain the values of {Os{^Si)) 
and {Os{^Po)) in the 2nd and 4th columns. The first errors in Table IV are extracted from 
the matrix of second derivatives of the function. There is also an error from varying the 
renormalization and factorization scales fip and fiR. These errors are large, but since they 
are highly correlated, we have separated them out as a second error in Table fV[ The upper 
and lower errors are the shifts in the matrix elements that minimize when fin and fip 



are changed from the central value /ly = + Pr by multiplicative factors of 2 and 1/2, 
respectively. The error from varying is also highly correlated, but it is smaller and we 
do not list it separately. It can be taken into account when we use the matrix elements to 
calculate other observables. 

Our fits for Bda / dpTdy at ?/ = for inclusive T(15'), T(2S'), and T(3S') are compared to 
the CDF data in Figures 1, 2, and 3. The error bands reflect the statistical uncertainties in 
the matrix elements. The fits are reasonably good in the region pt > 8 GeV that we used for 
fitting. At low pt, our fits for da/dpT behave like I/pt? because we have not implemented 
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3.6± i.9iJ:^ 

5.4±4.3l|^ 



3.9± i.7i?:^ 



5.7±4.6ty 


3.7± 2.lj:J;| 
7.5 ± 4.9t|^ 



4.1 ± i.9l:\i 


7.9 ± 5.2l|^ 




18.0 ± 5.611^ 
-10.2±9.7;?:8 



17.2 ± 5.0l|^ 


-10.6 ± 10.2;!^ 


19.6 ± 6.3t|^ 
-8.7 ± 11.1;? j 



19.0 ± 5.6l|J 


-8.9± 11.7;iJ 


b 


11.6 ±2.7+11 
10.9 ± 6.2^^°!^ 



12.4 ± 2.5t|^ 


11.1 ±6.5^105^ 


11.7±3.0l|:^ 
18.1 ±7.21^^^^ 



13.0 ± 2.8^1^ 


18.6 ± 7.5j:^^4^ 



TABLE IV. Inclusive color-octet matrix elements for T states (in units of 10 ^ GeV^). 

the effects of soft-gluon radiation on the shape of the pr-distribution. The fits therefore 
diverge from the data below = 8 GeV. For T{2S), the central curve becomes negative at 
small pt because our fit gives a negative central value for (O8(^5'o))j^c^'^'* or {08(^Po))JJc^^ ■ 
If we had fit the color-octet matrix elements using the data for px > 4: GeV, instead of only 
the data for p^- > 8 GeV, the central values for (O8(^5'o))inc or (08('^-Po))inc would also have 
been negative for T(IS') and T{3S). 

In Figures 1,2, and 3, the color-singlet model predictions from order-a^ fusion processes 
are shown as dashed lines. At the largest values of pt shown, the predictions fall more 
than an order of magnitude below the data. The color-singlet model prediction for T(3S') 
indicates the shape of the bbi{^Si) cross section. The color-singlet model predictions for 
T(liS') and T{2S) behave very differently at small pt, because they receive contributions 
from decays of Xbj{nP). The predictions diverge as px — >■ 0, because the cross sections 
da/dpT for &&i(^-Po,2) behave like 1/pt- In order to obtain the correct threshold behavior in 
these channels, it would be necessary to resum the effects of soft-gluon radiation. 



B. Direct Matrix Elements for P-waves 

The color-octet matrix elements (08('^>S'i)) for the Xfe's can be determined from CDF 
measurements of the fractions of T(lS')'s that come from Xf,(l-P)'s and from Xb(2P)'s [Q. 
The important feeddown decays for Xbj(l-P) and Xfej(2P) proceed through the T{2S) and 
T{3S), respectively. The fractions F^^^^p) of T(lS')'s from Xb{nP) therefore satisfy 

Fl^^% a[T(15)]i„, = E^x.,.(iP)^T(i5) cTlxuilP)] 

J 

+ (j2^'^i^s)^XbjiiP)Bxbj(iP)-'r{is)^ o-[T(2S')]i„c, (11) 

<(S cr[T{lS)U = E^x.,.(2P)-.T(i5) cr[xtj{2P)] 

J 
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FIG. 1. Inclusive cross section for T(IS') at y = in Run I multiplied by its branching fraction 
B into fj.~ as a function of pT- CDF data, NRQCD fit (solid line) with statistical error bars 
(dotted lines), and color-singlet model prediction (dashed line). 

+ (j2 ^'^i3S)-^xbj{2P)Bxij{2P)^r{is)^ a[T{3S)]inc- (12) 

Using the branching fractions from Table |, the coefficients of a[T{2S)]i^c and cr[T(3S')]inc 
in (|n|) and (0) are (3.8 ± 0.7)% and (1.8 ± 0.2)%, respectively. 

The CDF result for the fractions of T(lS')'s with \y\ < 0.4 and pt > 8 GeV that come 
from Xfe(lP)'s and Xb(2P)'s are 



^T{15) 
,T(15) 



(27.1 ±8.1) 



<(2P) = (10-5 ± 4.8) 



(13) 
(14) 



where we have added the statistical and systematic errors in quadrature. The inclusive 
T{nS) cross sections in ([TT|) and ([T2|) are the cross sections integrated over \y\ < 0.4 and 
Pt > 8 GeV. Using the CDF measurements of these cross sections, (|T3D and (0) reduce to 
the following constraints on the cross sections for Xbj{^P) and Xbj(2-P): 



J 



,P)^T{i5) ^[XbAnP)] = (0.85 ± 0.29) nb, n = 1, 
= (0.34 ±0.17) nb, n = 2. 
The branching fractions and the associated errors are given in Table |. 



(15) 
(16) 
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FIG. 2. Inclusive cross section for T(2S) at y = in Run I multiplied by its branching fraction 
B into fJ-^fJ-^ as a function of pT- CDF data, NRQCD fit, and color-singlet model prediction. 



The theoretical cross sections for a[xbj{nP)], J = 0,1,2, are obtained by integrating 
(0) over the appropriate region of y and Pt- The constraints (|T3p and ( ]TB| ) are then linear 
equations for (Of "^"^^(^^i)), which give the values in Table |V]. The first error is obtained 
by setting fip = f^R = f^T and combining in quadrature the experimental error from (^) or 
(P!B|) , the error from the branching fractions in Table |I|, and the errors from the color-singlet 
matrix elements in Table The second upper and lower errors in Table |V| are the shifts in 
the matrix elements when /ir and /ip are changed from their central values by multiplicative 
factors of 2 and 1/2, respectively. 



C. Direct Matrix Elements for S- waves 



The NRQCD matrix elements in Table ^ can be used to calculate the inclusive T{nS) 
cross sections. To calculate the direct T{nS) cross sections, we must extract direct color- 
octet matrix elements for the T{nS) states from the inclusive color-octet matrix elements 



given in Table IV. The linear combinations of matrix elements determined by the inclusive 



T(IS') cross sections are 

{OsCSi))l^^'^ = {Oj^''\'S,)) + 0.311 {Oj^''\' Si)) +0.112 {Oj^''\'Si)) 

+2.15 (Or^'^^H'^i)) + 1-08 (Or^'^'^e^i)), (17) 
{OsCSo))!^''^ = {Oj^''\'So)) + 0.311 {Oj^''\'So)) + 0.112 {Oj^''\'So)), (18) 
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{OsCPo))l^'^ = {Oj^''\'Po)) + 0.311 {Oj^''\'Po)) + 0.112 {0^^''\'Po)). (19) 

The errors on the branching fractions in (|T7|)-([l9|) have been suppressed, but they are 
(0.311 ± 0.016), (0.112 ± 0.005), (2.15 ± 0.31), and (1.08 ± 0.08). The hnear combinations 
determined by the inclusive T(2S') cross sections are 

{Os{'S^)JS'^ = {0^^''^eS,)) + 0.106 (Oj('^)(35i)) + 1.49 {0^^'''^' S,)) , (20) 
{OsCSo))JS'^ = {Oj^''\'So)) + 0.106 (Oj^^^V^o)), (21) 
{OsCPo))IS'^ = {Oj^''\'Po)) + 0.106 {0^^''\'Po)). (22) 



The errors on the branching fractions in (|0D-(|2D are (0.106 ± 0.008) and (1.49 ± 0.17). 
Using the color-octet matrix elements for inclusive T[nS) in Table |V| and the color-octet 
matrix elements for direct Xb{nP) in Table 0, we obtain the color-octet matrix elements for 
direct T{nS) in Table |V|. 

Our analysis gives a negative value consistent with zero for the matrix elements 
{Os^^^\'So)) or {Oj^^^\''Po)). Our values for (©^^'''^(^Si)), {0^''^^''\^Si)), and 
{Og^^^\'^Si)) are also consistent with zero given the statistical error. The only direct 
color-octet matrix elements that differ from zero by two or more statistical error bars are 
{Oj^''\'S,)) and {Oj^''\'So)) or {0^^''\'Po)). 

We now compare our values for the matrix elements with those obtained by Cho and 
Leibovich in their pioneering analysis of bottomonium production at the Tevatron. Their 
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CTEQ5L 


MRSTLO 






i-2±i.3tJ:^ 




1.5 ±1.111:3 


i.9± i.3tl:^ 




16.4±5.7t^:} 
-10.8±9.7;3:^ 



15.6 ± 5.21^:^1 


-11.2 ± 10.2;!^ 


17.4 ±6.4;^:? 
-9.5±ll.i;|8 




16.8 ± 5.8tti 


-9.7±11.6;|? 


b 


2.0±4.1+[J:6 
13.6 ± 6.8tl%^ 



3.0 ± 3.8j:[j:? 



13.9±7.lt^^o^ 


0.4 ±4.7;!:? 

20.2 ± 7.81^^5^ 



1.8±4.4;[J:2 


20.7 ±8.1+^1^ 



TABLE V. Direct color-octet matrix elements for Xb ^-nd T states (in units of 10 ^ GeV^). 



analysis was based on the CDF data sample from Run lA of the Tevatron [§. To reduce 
the errors associated with the shape of the pt distribution at small pt, they used only the 
data from p^ > 3.5 GeV in their analysis. The data was insufficient to determine all the 
matrix elements, so they estimated the matrix elements {0^^"'^\^Si)) by applying scaling 
relations to the corresponding matrix elements in the charmonium sector. Their value for 
{Oj^^^\^Si)) is consistent with ours to within our large error bars, but their value for 
{Og^'^^\^Si)) is smaller than ours by about a factor of 40. They used the CDF data to fit 
the matrix elements {Og''^"^\^Si)) and the linear combinations M^^"'^^ = {Oj^"'^\^So)) + 
5{0g {^Po)) /ml. Their values for {Og i^^i)) are comparable to ours in magnitude, 
but they have much smaller error bars. Their values for mJ^""^^ differ from zero by only 
about one error bar. In our analysis, we included the matrix elements {Oj''"'^\^Si)) in the 
list of those to be fit to the CDF data. The much higher quality of the CDF data from Run 
IB of the Tevatron 0] allowed us to carry out a reasonable fit using the data restricted to 
Pt>8 GeV. 



VI. INCLUSIVE CROSS SECTIONS FOR SPIN-TRIPLET STATES 

Having determined the most important matrix elements for the production of the spin- 
triplet bottomonium states, we can use them to calculate the cross sections for these states 
in other high energy processes. In particular, we can calculate their cross sections in Run 
II of the Tevatron in which the center-of-mass energy will be increased from 1.8 TeV to 2.0 
TeV. In order to cancel the large theoretical errors, such as those from the uncertainties in 
the matrix elements and from the choice of scale, we normalize the cross sections to that for 
inclusive T(IS') at ^/s = 1.8 TeV. For each bottomonium state H, we define the ratio 

^""^^^ " ^[inclusive T(l^);^ = 1.8 TeV] ' ^^^^ 
where the cross sections are integrated over pt > 8 GeV and over \y\ < 0.4. 
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H 


R^{1.8 TeV) 


R"{2.0 TeV) 


T{3S) 


0.31 ±0.14 


0.36 ±0.16 




0.44 ± 0.26 


0.52 ±0.30 


Xbi{2P) 


0.34 ±0.16 


0.39 ±0.19 




0.20 ± 0.07 


0.24 ± 0.08 


T{2S) 


0.65 ±0.35 


0.76 ±0.41 


Xb2{lP) 


0.57 ±0.26 


0.66 ±0.31 


XbiilP) 


0.41 ±0.17 


0.48 ±0.19 


XboilP) 


0.23 ± 0.08 


0.26 ±0.09 


T{1S) 


1 


1.16 ±0.01 



TABLE VI. Ratios of the inclusive cross sections for the spin-triplet bottomonium states H at 
the Tevatron with ^/s = 1.8 TeV and 2.0 TeV to the inclusive cross section for T(IS') with ^/s = 
1.8 TeV. 



To calculate the inclusive cross section for T{nS), we simply use the inclusive color-singlet 
matrix elements from Table |T| and the inclusive color-octet matrix elements from Table 
To calculate the inclusive cross section for Xbj{nP), we must first compute the direct cross 
sections for Xbj{nP) and the higher bottomonium states using the direct color-singlet matrix 
elements from Table and the direct color-octet matrix elements from Table 0, and then 
combine them using the inclusive branching fractions in Table |. The resulting ratios 
shown in Table [V^ are the averages of the 4 values obtained by using either the CTEQ5L 
or MRST98LO parton distributions and either setting (OsCSq)) = or {Os{^Po)) = 0. The 
error bars come from combining in quadrature the statistical errors in the matrix elements 
from the Tables, the error from varying yU by a factor of two from its central value, the 
difference between using the CTEQ5L and MRST98LO parton distributions, the difference 
between setting (Osi^So)) = and (Osi^Po)) = 0, and the error from varying rrib. The 
error bars in the numerator and denominator of are both large, but they are highly 
correlated and tend to cancel in the ratio. The largest contributions to the error bars are the 
statistical errors on the matrix elements, with the exception of T(IS'), for which the largest 
contribution comes from varying the scale. For T{2S) and T{3S), the results in Table ^ 
for ^/s = 1.8 TeV are consistent with the actual CDF measurements, which give 0.61 ±0.12 
for T{2S) and 0.29 ± 0.12 for T(3S'). When the center-of-mass energy is increased from 1.8 
TeV to 2.0 TeV, all the cross sections increase by about 16%. The increase depends on pr, 
changing from about 15% at pr = 8 TeV to about 19% at pr = 20 TeV. 



VII. DIRECT CROSS SECTIONS FOR SPIN-SINGLET STATES 



Having determined the most important matrix elements for the production of the spin- 
triplet bottomonium states, we can also use them to calculate the production rate of the 
spin-singlet states rif,{nS) and hh{nP). The matrix elements for these states are related to 
those of the corresponding spin-triplet states by the approximate spin symmetry of NRQCD. 
Spin symmetry relates the matrix elements for rib{nS) to those for T{nS): 
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T(n5)/3 



Si)), 



{Or^'\'Po)), 
{Ol^'''\'S,)), 



;-{Oj^'''\'So)). 



(24) 
(25) 
(26) 
(27) 



The direct matrix elements for rii,{nS) can therefore be read off from those for T{nS) in 
Tables ^ and |V|. The spin-symmetry relations have been used to calculate the cross sections 
for producing the r/c at the Tevatron [|16|, in photoproduction and electroproduction, [ll7| , 
and at Hera-B [|18|]. In Refs. and |T^, the factor of | which comes from the number of 
spin states was omitted from the spin-symmetry relations (p^)- (|27|) . Thus their estimates 
of the cross sections may be too large by a factor of 3. Spin symmetry relates the matrix 
elements for hb{nP) to those for Xbo{nP): 



{Ol''^-P\^P^)) = 3{0^^o{n''\^Po)), 



(28) 
(29) 



The direct matrix elements for hf,{nP) can therefore be read off from those for Xb{nP) in 
Tables ^ and |V|. These relations were first used by Fleming and Mehen to calculate the 
cross section for the photoproduction of he ||19|- They have also been used to calculate the 
cross sections for he at Hera-B ||18|| . 

To calculate the cross sections for the direct production of the spin-singlet states at the 
Tevatron, we need the appropriate parton differential cross sections da. Explicit expressions 
for most of those that are needed are given in Ref. [Q] and in Ref. 0. The exception is the 
cross sections for producing 663 -Pi), which contributes to rjf, production. The cross sections 
for the production of bbs{^Pi) in qq, gq, and gg scattering are 



da 
di 
da 
di 
da 



{qq ^ 'P^g 
{gq ^ g 



{Osi'Pi))AiAna,)^ ? + 

167rs2 9M3 s(s-M^)^'' 
(Osi'Pi)) (47r«,)3 



IQttP QM^ (_t)(M2-t)2' 

{Osi'Pi)) (47ra,)3 1 

167rs2 36M3 z^rs{rs - M^)^{sM^ + z^) 



(30) 
(31) 



{275^2(58 - 4sV + + + ^«) 

+M2(27s^° - 2A3s^z^ + Q97s'^z^ - 665s^z^ + 3465^/ - 27^^°) 
-M^s(135s^ - 7025^^2 + 1340sV - 10875^2^ + 135z^) 
+M^(324s* - llSis^z^ + 1557 s^z^ - 698s V + 54z^) 
-M*s(486s^ - 1091sV + 8825^^^ - 92z*^) 
+M^°(486s^ - 616sV + 3745^/ - 27/) 



M^^s(324r - 211s^z^ + 38z 



+M^\sHl35s'^ 



38z' 



27M 



16-3 



(32) 
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±1 


It \^i.o le V ) 


it le V ) 


%(35) 


1.72 ±0.52 


2.00 ±0.61 


hb{2P) 


0.07 ±0.07 


0.08 ± 0.09 




1.77 ±0.50 


2.06 ±0.59 


hilP) 


0.11 ±0.08 


0.13 ±0.09 


VbilS) 


4.31 ±0.98 


5.02 ± 1.14 



TABLE VII. Ratios of the direct cross sections for the spin-singlet bottomonium states H at 
the Tevatron with ^/s = 1.8 TeV and 2.0 TeV to the inclusive cross section for T(15) with y/s = 
1.8 TeV. 



where =iu. 

We proceed to calculate the cross sections for the direct production of the spin-singlet 
states at the Tevatron at center-of-mass energies = 1.8 TeV and 2.0 TeV. To minimize 
the effect of the highly correlated errors, we calculate the ratio (p3D of the direct cross 
section integrated over > 8 GeV and over \y\ < 0.4 to the corresponding inclusive cross 
section for T(IS') at 1.8 TeV. The resulting predictions are shown in Table 0. The cross 
sections for the hf,{nP) states are small compared to those for T(IS') and they have large 
error bars. The cross sections for r]h{nS) are predicted to be several times larger than those 
for T{nS) and they have reasonably small error bars. These predictions should be fairly 
reliable, because the cross sections are for the same kinematical region as the data used to 
extract the matrix elements. When the center-of-mass energy is increased from 1.8 TeV to 
2 TeV, all the cross sections increase by about 16%. 

We can make a rough estimate of the cross sections integrated over all px by assuming 
that the spin-singlet cross sections have the same shape at small pr as the T(IS') cross 
section. The measured inclusive T(IS') cross section at central rapidity integrated over 
all pt up to 20 GeV satisfies Bda/dy = 690 ± 25 pb, where B ^ 2.5% is the branching 
fraction of T(IS') into fi'^fi~. The cross section integrated only over p^ > 8 GeV satisfies 
Bda/dy = 106±7 pb. The ratio of these cross sections is 6.5±0.5. Multiplying the inclusive 
T(IS') cross section da/dy = 28 nb by the factor of 6.5 to take into account the small pt 
region and by the ratio 4.3 from Table we find that the cross section for rii,{lS) integrated 
over all pt should be approximately da/dy = 800 nb. 

The cross section for rjbi^S) indicates that this state must have been produced in abun- 
dance in Run I of the Tevatron. However the rii,{lS) can be observed only if it has a large 
enough branching fraction into a decay mode that can be triggered upon. One possibility is 
the double- J/?/' decay r]h{lS) J/ip + J/'ip, followed by the decays J/ip ^ fi~^fi~ . The decay 
rib{lS) — >• J /ip + J /ip has essentially the same kinematics as the decay rjc — > 00, except that 
all masses are scaled up by a factor of 3. Thus the branching fraction for rib{lS) J /ip+J/ip 
could be as large as that for rjc 00, which is approximately 7x 10^^. We can obtain a lower 
bound on the branching fraction by using the fact that in the limit rrib oo with rric fixed, 
the branching fraction for rjb{lS) — > J/ip + J/ip scales like 1/mf ||2^. If rric and rrib were both 
in this scaling region, then the branching fraction into light J^'^ = 1 mesons would be 
smaller for rjbi^S) than for r/c by a factor of {rric/mbY, which is about 10"^. Since we are not 
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deep into this scaling region, the suppression should be smaller than this. Thus the branching 
fraction for ?7fe(lS') — >■ J /ip + J /tl) should be in the range between 7 x 10"^ and 7 x 10"^ Mul- 
tiplying by the branching fractions of 6% for each of the decays J /ip ^ f^~^f^~, our estimate 
for the branching fraction for ?7f,( IS") — > J/ip + J/ip — fi~^ fi~ + fi~^ fi~ is i? ~ 2.5x10"^^^. The 
cross section for producing this particular decay mode of rif){lS) is therefore Bda/dy ~ 2 pb, 
give or take a factor of 10. Multiplying by the rapidity interval 0.8 and by the integrated 
luminosity of about 100 pb^^ in Run I of the Tevatron, we obtain between 16 and 1600 
produced events. 

We must also take into account the acceptances and efficiencies for observing the decays 
J/ip ^ fj'~^fi~- These can be estimated using the CDF data on the production of prompt 
J/ip in Run lA of the Tevatron Based on the observation of about 22,000 J/ip —* ^^fi~ 
candidates with pt > 5 GeV and pseudorapidity |?7| < 0.6 in a 15 pb~^ data sample, they 
measured the cross section in that region of pt and t] to be Ba ~ 17 nb. We infer that the 
product of the acceptance and the efficiency is roughly e ~ 0.09. Multiplying the number 
of events that are produced by e^, we get between 0.13 and 13 observed events. Thus 
this cross section may be large enough to be observed in Run I of the Tevatron. In Run 
II, the integrated luminosity will be larger by a factor of 20 and there will be significant 
improvements in the acceptances and efficiencies for observing muons in both the CDF and 
DO detectors. Thus the ?7b(15') should certainly be observed in Run II through the decay 



VIII. DISCUSSION 

We have carried out an updated NRQCD analysis of the CDF data on the production 
of spin-triplet bottomonium states from Run I of the Tevatron. In spite of using only the 
data from pt > 8 GeV, we were able to extract all the relevant color-octet matrix elements 
directly from the data. Only one of the 8 color-octet matrix elements comes out with a 
negative central value, but several others are also consistent with zero to within errors. In 
our analysis, we distinguished between the inclusive color-octet matrix elements that can be 
used to compute inclusive T{nS) cross sections and the direct color-octet matrix elements 
required to compute direct T{nS) cross sections and, by spin symmetry, direct r]f,{nS) cross 
sections. 

The most serious deficiency in our analysis was our failure to take into account the effects 
of soft-gluon radiation that are needed to give a smooth pj- distribution at small pt- This 
forced us to use only the small fraction of the data from p-r > 8 GeV to ffi the color- 
octet matrix elements, which led to large errors in these matrix elements. If these matrix 
elements are used to predict bottomonium cross sections in other high energy processes, the 
predictions will probably have large error bars. Our theoretical cross sections also diverge 
from the CDF data below pt = 8 GeV, which gives us another reason to be cautious in 
applying our matrix elements to other high energy processes. 

An analysis that deals properly with the small pt region could take full advantage of the 
CDF data and therefore determine the color-octet matrix elements more accurately. Such an 
analysis requires a prescription for combining the leading-order cross sections for ij —>■ bb + k 
with the next-to-leading order cross sections for ij bb recently calculated by Petrelli et al. 
T2\ to get a smooth px distribution near px = 0. The matrix elements extracted from such an 
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analysis should give reliable predictions for observables at the Tevatron that are dominated 
by low pt- They should also give reliable predictions for bottomonium production in other 
high energy processes. 

Our analysis should give reliable predictions for the cross sections of the spin-singlet 
states r]b{nS) and hb{nP) at the Tevatron. We find that the direct cross section for rjb{lS) 
at Pt > 8 GeV should be greater than the inclusive cross section for T(IS') by a factor of 
about 4.3. If we assume that the pt distributions have the same shapes at smaller pt, we 
can estimate the cross section for rjbi^S) integrated over all px- The resulting cross section 
is large enough that it should be possible to discover the r]b{lS) at the Tevatron if there 
is a decay mode with a large enough branching fraction that can be used as a trigger. We 
argued that the decay rjb J/ip + J/ip should allow the discovery of the r]b{lS), if not in 
the data from Run I, then certainly in Run II of the Tevatron. 
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